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We propose a global measure for quantum correlations in multipartite systems, which is obtained 
by suitably recasting the quantum discord in terms of relative entropy and local von Neumann 
measurements. The measure is symmetric with respect to subsystem exchange and is shown to be 
non-negative for an arbitrary state. As an illustration, we consider tripartite correlations in the 
Werner-GHZ state and multipartite correlations at quantum criticality. In particular, in contrast 
with the pairwise quantum discord, we show that the global quantum discord is able to characterize 
the infinite-order quantum phase transition in the Ashkin- Teller spin chain. 



PACS numbers: 03.67.-a, 75.10.Pq, 03.67.Mn 



I. INTRODUCTION 

Quantum correlations constitute a fundamental re- 
source for quantum information tasks pj. They are 
rooted in the superposition principle, displaying effects 
with no classical analog. The research on quantum cor- 
relation measures was initially developed based on the 
entanglement-separability paradigm [2]. More recently, 
however, it has been perceived that entangled states are 
not the only kind of quantum states exhibiting non- 
classical features. In this context, a suitable measure 
of quantum correlation has been introduced by Ollivier 
and Zurek Q. This measure, which has been designated 
as quantum discord, is able to capture not only quan- 
tum correlations in entangled states but also in separable 
states. It arises as a difference between two expressions 
for the total correlation in a bipartite system (as mea- 
sured by the mutual information), which are classically 
equivalent but distinct in the quantum regime. Remark- 
ably, quantum discord has been revealed as a useful quan- 
tity in a number of applications, such as quantum crit- 
ical phenomena [U, H| and quantum evolution under de- 
coherence [6]. Moreover, quantum discord has also been 
conjectured to be a resource for speed up in quantum 
computation 0] and for locking classical correlations in 
quantum states @. 

In recent years, generalizations of quantum discord to 
multipartite states have been considered in different sce- 
narios ||. One possible approach is based on directly 
generalizing the quantum mutual information to a mul- 
tipartite system, even though nonunique generalizations 
are possible in this situation [l(J [H| • Another approach 
is to define from the beginning a measure based on the 
relative entropy, which allows for a unified view of differ- 
ent correlation sources, such as entanglement, quantum 
discord, and dissonance [IH (see also Ref. [Hj). The 
aim of this work is to propose a global measure of quan- 
tum discord obtained by a systematic extension of the 



bipartite quantum discord, with operational appeal and 
satisfying the basic requirements of a correlation func- 
tion. In this direction, we suitably recast the standard 
bipartite quantum discord defined in Ref. [3| in terms of 
relative entropy and local von Neumann measurements, 
whence a natural multipartite measure for quantum cor- 
relations emerges. This measure - named here as global 
quantum discord (GQD) - is symmetric with respect to 
subsystem exchange and shown to be non-negative for ar- 
bitrary states. We illustrate our results by computing the 
tripartite GQD in the Werner-GHZ state and by applying 
GQD in the characterization of the infinite-order quan- 
tum phase transition (QPT) in the Ashkin- Teller spin 
chain, where ordinary pairwise quantum discord fails. 



II. QUANTUM DISCORD 

Consider a bipartite system AB composed of subsys- 
tems A and B. Denoting by pab the density operator of 
AB and by pa and ps the density operator of parts A 
and B, respectively, the total correlation between A and 
B is measured by the quantum mutual information 



I (pab) = S(p A ) - S(pa\pb), 



(1) 



where S(pa) — —TrpA l°g 2 Pa is the von Neumann en- 
tropy for A and 



S(pa\pb) = S(pab) - S(p B ) 



(2) 



is the entropy of A conditional on B. The conditional 
entropy can also be introduced by a measurement-based 
approach. Indeed, consider a measurement locally per- 
formed on B, which can be described by a set of pro- 
jectors {I1 B } = The state of the quantum 
system, conditioned on the measurement of the outcome 
labeled by j, is 



PAB\ 
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Pab 



(3) 
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where pj = Tt[(1a <8> H 3 b )pab(1a ® Hg)] denotes the 
probability of obtaining the outcome j and 1a denotes 
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the identity operator for A. The conditional density op- 
erator Pab\j allows for the following alternative definition 
of the conditional entropy: 

S(p AB \{I! j B }) = YlPi S (P^ W 

3 

where p A \j = Tt b PAB\j = 0-/Pj)(bj\pA B \bj), with 
S(pA\j) = S(pAB\j)- Therefore, the quantum mutual in- 
formation can also be defined by 

J(p AB ) = S(p A ) - S(p AB \{& B }). (5) 

The quantities I (pab) and J(pab) are classically equiv- 
alent but they are distinct in the quantum case. This 
difference is the quantum discord T> (pab) [3j, yielding 

T>{pab) = I(pab) - J (pab) ■ (6) 

Note that T> (pab) is defined as a non- negative asymmet- 
ric quantity that depends on {II S }. This dependence can 
be eliminated by minimizing T> (pab) over all measure- 
ment bases {fl j B } [3]. 

III. RELATIVE ENTROPY AND SYMMETRIC 
QUANTUM DISCORD 

The quantum relative entropy is a measure of distin- 
guishability between two arbitrary density operators p 
and a, which is defined as j 1 51 ] 

S(p || a) = Tr(plog 2 p-plog 2 <7) . (7) 

We can express the quantum mutual information I(pab) 
as the relative entropy between pab and the product 
state pa® Pb, i-e. 

I (pab) — S (pab II Pa €D pb) ■ (8) 

In order to express the measurement-induced quantum 
mutual information J (pab) in terms of relative entropy, 
we need to consider a non-selective von Neumann mea- 
surement on part B of pab , which yields 

$b (pab) = ® ^b) Pab (Ia ® n J B ) 

3 

= J2P3PA\j®\b j ){b j \. (9) 

3 

Moreover, tracing over the variables of the subsystem A, 
we obtain 

$s (Pb) = $b (Tr A pab) = \bj) (b s \ , (10) 

3 

where we have used that Tr A (p A u) — 1. Then, by ex- 
pressing the entropies S (&b (pab)) and S ($ B (pb)) as 

S ($b (pab)) = H (p) + J^PjS (pa\j) (11) 

i 



and 

S($ B (p B )) = H(p), (12) 
with H (p) denoting the Shannon entropy 

ff(p) = -5>ilog 2 ( Pj ), (13) 

3 

we can rewrite J (pab) as 

J (pab) = S (p A ) - ^PjS (p A \j) 
j 

= S (p A ) + S ($ s ( PB )) - S (<S> B (pab)) 

= S($b(pab) || pa®^b(pb))- (14) 

Therefore, the quantum discord can be rewriten in terms 
of a difference of relative entropies: 

£> (pab) = S (pab II PA <8> Pb) 

-S($ b (pab)\\ Pa®$b(pb)), (15) 

with minimization taken over {II B } to remove the 
measurement-basis dependence. It is possible then to 
obtain a natural symmetric extension V (pab) for the 
quantum discord V(pab)- Indeed, performing measure- 
ments over both subsystems A and B, we define 

D (pab) = min [S (pab \\ PA <S> Pb) 
{n^®n*} 

-S ($ab (pab) II ®A (pa) <8 $b (ps))] , (16) 
where the operator <I>ab is given by 

$AB (PAB) = { & A ® PAB (ft J A ® ft|) • (17) 

J,* 

Observe that, by writing Eq. (|16p in terms of the mutual 
information /, we obtain 

V (pab) = min [I {pab) - I($ab (Pab))] , (18) 

which is the symmetric version of the expression for the 
loss of correlation due to measurement [ll|, [lj| . Remark- 
ably, 2? (pab) is equivalent to the measurement-induced 
disturbance (MID) [I?) if measurement is performed in 
the eigenprojectors of the reduced density operators of 
each part (instead of minimization). Moreover, Eq. (|16p 
also provides the symmetric quantum discord considered 
in Ref . [l8[ and experimentally witnessed in Ref . [ll| . As 
a further step, we can still rearrange Eq. (|16[) in a rather 
convenient way, yielding 

£> (pab) = min [S (pab || $ab (pab)) 
{n J A ®n k B } 

-S (p A II $a (pa)) - S (p B I) $ B (p fl ))] • (19) 
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IV. GLOBAL QUANTUM DISCORD 

Let us now extend quantum discord as given by 
Eq. (fP9]) to multipartite systems. 

Definition. The global quantum discordV (p Al ... AN ) for 
an arbitrary multipartite state p Al ... AN under a set of 
local measurements {H Al ® ■ ■ ■ <g> } is defined as 

T>(p Al -A N ) = min[S (p Al ---A N || $ (p Al ... Aff )) 
{rM 

N 

-J2 s (p*i IIMp^))L (20) 

3=1 

where <f>j (p Aj ) = £\, n A, PA, tl A . and $ ( PAi ... An ) = 

Y, k ftfc PA X -A N flfc, with Il k = fl A <8 • • • ® ft^. a«rf k 
denoting the index string {j\ • • ■ Jn )■ 

Therefore, a classical state can be denned by 
PAi---A N — & {pAx-'-An), which is in agreement with 
the requirement that classical states are not disturbed 
by suitable local measurements. Indeed, this defini- 
tion of a classical state implies that pa, = &j (pa,) for 
any j, which means V (p Al ... AN ) — 0. Moreover, ob- 
serve that, via minimization over the set of projectors 
{il Al ® • •• ® n-^}, we define GQD as a measurement- 
basis independent quantity. However, as will be illus- 
trated in the Ashkin- Teller chain, other (non-minimizing) 
bases are also able to provide relevant information about 
the behavior of quantum correlations in the system (sim- 
ilarly to the original definition of quantum discord in 
Ref. [3j). In any case, we can show that GQD is non- 
negative for an arbitrary state. 

Theorem. The global quantum discord T> (p Al ... AN ) is 
non-negative, i.e., T> (p Al ... AN ) ^ 0. 

Proof. In order to prove that T> (p Al ... AN ) ^ 0, we as- 
sociate with each subsystem A, an ancilla system Bj. 
Therefore, we will define a composite density operator 

PAx-AmBi-Bn such that 

p'\ —A \ :/; = 5353 ^-kPAt-A N ftfc' ®hk>, (21) 
k k' 

where k kk , = \B n )(B fi \<E) •••(g) \B jN )(B f J, with k and 
k' denoting the index strings (ji • • ■ Jn) and (j[ ■ ■ ■ j' N ), 
respectively. From the monotonicity of the relative en- 
tropy under partial trace (2(j, for any positive opera- 
tors <Ji2 and 712 such that Tr(<7i2) = Tr (712)1 we have 
that S , (<5-i 2 ||7i2) > 5((Ti||7i), where a\ = Tr 2 (<7i 2 ) 
and 7i = Tr 2 (7i 2 ). Then S (ai23...Ar||7i23..Jv) > ■ •■ > 
S (0123II7123) > -5" (0-12II712) ^ 5(o-i||7i). By taking 

P , A 1 --A N ;B 1 -B N aS ° and PA X ;B X ® p'a 2 ,b 2 ® • • • ® Pa n ;B n 

as 7, we obtain 

S {PA-l-An^-B* WP'auBi ® P'a 2 ;B 2 ® ■ ■ ■ ® PAjviBiv) 

>S'(/3a 1 ...a w IIA 1 ®--®P^). (22) 



which, from Eq. ((SJ, implies that 

iV 

(PA,;B,) ~ 5 (PAi-AjviBi- -Sjv) 

3=1 

JV 

>12 S {pa 3 )- S (p'a 1 ...a n )- (23) 
3=1 

However, from Eq. (f2T|) . it follows the relations 

5 (PAi"' J 4jv;Bi"'Bjv) = S {pA t -A N ) , (24) 

S(p'a 1 -a n )=S(<Z>(pa 1 -a n )) , (25) 
S (^., Bj ) = S (p Aj ) (Vi), (26) 
=5(^(^)) (Vj). (27) 

Insertion of Eds. ([21 ]) - ([2?] ) into inequality yields 

JV 

53 5 (^i) ~ S (PAi-A N ) 
3 = 1 

AT 

^J2 S ^ApA j ))-S(^(p Al ... AN )) . (28) 
3+1 

By rewriting inequality l[28p in terms of the relative en- 
tropy, we obtain 

AT 

S(fi Al ...AA^(fiA 1 ...A ll ))-^S(pA i \\ij{pA i ))>0. 

3=1 

(29) 

The left hand side of the inequality above is exactly the 
GQD, as defined by Eq.®. Hence, V {p Al -A N ) > 0. ■ 



V. TRIPARTITE CORRELATIONS IN THE 
WERNER-GHZ STATE 

As a first illustration of GQD, we will consider the 
Werner-GHZ state 

p = C 1 ~ P) l + ^\qhz) (GHZ\ , (30) 
8 

where < // < 1 and 

\GHZ) = {\\) A \\) B |t) c + ||> , 4 ||) B \i) a )/V2, (31) 

with |t) and \\.) denoting the eigenstates of the Pauli 
operator a z associated with eigenvalues 1 and — 1, re- 
spectively. The Werner-GHZ state provides an interpo- 
lation between a fully mixed (uncorrelated) state and a 
maximally correlated pure tripartite state. It is a rather 
suitable state to begin with as we propose a measure 
for quantum correlation and constitutes an interesting 
scenario to compare multipartite with bipartite correla- 
tions, since it is a generalization of the two-qubit Werner 
state [2l|. Let us begin by analyzing GQD in the case of 
a pure GHZ state (p = 1). 
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A. GQD for the GHZ state 

Let us focus here on the GHZ state, as defined by 
Eq. (|3"Tj) . In order to define local measurements for 
\GHZ), let us consider rotations in the directions of the 
basis vectors of subsystems A, B, and C, which are de- 
noted by 



We plot D(#2,#3) as a function of # 2 and #3 in Fig. Q] 
Notice that its minimum, which provides T> (p) , occurs 
at the boundary values 82 = #3 = 0, where T> (p) — 
1. This is a manisfestation of the fact that any local 
measurement disturbs the GHZ state, which is detected 
by a nonvanishing GQD. 



t) 3 



(32) 



|), , (33) 



with j = 1,2,3 for subsystems A, B, and C, respectively. 
The angles 9i take values in the interval [0, ir) and the 
angles ifi take values in the interval [0, 2n). In order to 
compute T> (p), with p — \GHZ) (GHZ\, we must evalu- 
ate the expression 



V (p) = min [S(p\\ <&(p)) 
-S (pi 



S (pa II $4 (pa)) 
$b (p B ))~S(p c II *c(pc))] -(34) 



However, S (p) = 0, since p is pure. More- 
over, S (pA II ®A (PA)) = S (p B || ®B (PB)) 

S (pc II ®c (pc)) = 0, since p^, p B , and p c are 

proportional to identity operators. Hence, GQD is 
simply given by 



V(p) 



min S ($ (p)) = min 



j log 2 



(35) 



CO 

CD 
<N 
O 




1 9 3 (rad) 



9 2 (rad) 



FIG. 1. The function 0(62,63) as a function of 62 and S3. 
Note that the minimum occurs for 62 = 63 — 0, which implies 
V(p) = l. 



where Xj are the eigenvalues of the operator $ (p) . They 
can be obtained from projections of the GHZ state over 
the rotated basis states. In order to minimize S (<& (p)), 
we must find out the measurement basis that maximizes 
the purity of $ (p), i.e., that maximizes the dispersion of 
the eigenvalues Xj with respect to the average of {Xj}. 
This is obtained for Bi — (i=l,2,3), namely, measure- 
ments in the eigenprojectors of cr? . As an illustration, let 
us consider the case of Q\ = and <p.; = (i=l,2,3). In 
this situation, the eigenvalues Xj for the operator $ (p) 
read 



Ai — As — 
A2 = A7 = 
A3 = A6 = 
A4 = Ah 



- cos — COS . 
2 V 2 J V 2 



-cos \- sm ( y 



1 ■ 2 / ^2 \ 2 / ^3 

-sm - cos U- 



1 ■ 2 / p 2 \ . 2 / ^3 

- sin — sin — 

2 \2 2 



(36) 
(37) 
(38) 
(39) 



By using Eqs. (H)])-® into Eq. ([35]), we can directly 
obtain T> (p) by minimizing over the angles 62 and #3. 
The function D(02, 63) to be minimized is then 



D(6o,e a 



zZ X 3 lo S 



2 Aj- 



(40) 



B. GQD in the Werner-GHZ state 

In order to obtain V (p) for the Werner-GHZ state, let 
us first rewrite Eq. (|3"0|) as 



P 



'1"3 



'2 U 3 



'1"2"3 



-ololol 



<J 1 (T 2 <T 3 



(41) 



with erf, o\ , and of denoting the Pauli matrices for the 
qubit i. Again, we will have here that S (pi || <!>i (pi)) = 
(i = A, B, C), since pa, Pb, and pc are proportional to 
identity operators. Therefore, 



V (p) = min 5 (pp (p)) = min [5 ($ (/>)) - 5 (p)] 

The von Neumann entropy S (p) is given by 
7 



(42) 



S(p)=3--(l-/i)log 2 (l- M ) 



1 



(l + 7p)log 2 (l + 7p) 



(43) 



As for the GHZ state, we take local measurements in 
the <7 Z eigenbasis for each particle to minimize S ($ (p)). 
Such an eigenbasis provides the maximum loss of corre- 
lation among the parts of p, which therefore minimizes 



5 



GQD. Then, from Eq. ([3Pp. we obtain 



1 — /i \ - /i 



8 



which implies 



(44) 



S($(p)) = 3- 3 (l- / x)log 2 (l-/x) 

-i(l + 3/x)log 2 (l + 3/i) . (45) 
Insertion of Eqs. (g3]) and (gS|) into Eq. (g2]) yields 
P(/5) = -i(l + 3/i) log 2 (l + 3 M ) 

+ i (1 - M ) log 2 (1 - ^ + 1 (1 + 7 M ) log 2 (1 + 7 M ) .(46) 

In Fig. [2] we plot 2? (p) as a function of /i. Observe that 
GQD vanishes only for p = 0, where p is a completely 
mixed state. Moreover, GQD is a monotonic function of 
p, acquiring its maximal value T> (p) — 1 for p = 1, where 
p is the GHZ state. This result resembles the behavior 
of the bipartite Werner state [3, H3| • 




FIG. 2. (Color online) Tripartite GQD for the Werner-GHZ 
state as a function of the mixing parameter p. Note that 
GQD is nonvanishing for p ^ 0. 



VI. MULTIPARTITE CORRELATIONS IN THE 
ASHKIN-TELLER CHAIN 

Let us now present an application of GQD that makes 
evident the importance of considering genuine multipar- 
tite correlations to the characterization of a QPT. In this 
direction, we consider the Ashkin- Teller model, which 
has been introduced as a generalization of the Ising spin- 
1/2 model to investigate the statistics of lattices with 
four-state interacting sites |23j . It exhibits a rich phase 
diagram 24j and has recently attracted a great deal 
of attention due to several interesting applications [25$ . 



The Hamiltonian for the quantum Ashkin- Teller model 
in one-dimension for a chain with M sites is given by 



M 



M 



-^/^E(^ + i+^^i+ A ^^ + i^i)- (47) 

3=1 

where a" and f? (a — x, y, z) are independent Pauli 
spin-1/2 operators, J is the exchange coupling constant, 
A and /3 are (dimensionless) parameters, and periodic 
boundary conditions (PBC) are adopted, i.e., <7^ +1 = <r" 
and Tm+i = 'Tt* ( a ~ x tV-, z )- The Ashkin- Teller model 
is Z2 ® Z2 symmetric, with the Hamiltonian commuting 
with the parity operators 



Pi 



and 



V 2 



(48) 



Therefore, the eigenspace of Hat can be decomposed 
into four disjoint sectors labeled by the eigenvalues of V\ 
and V2, namely, Q = {V x = +1,V 2 = +1), Q = 1 
{Vx = +1,V 2 = -1), Q = 2 {Vi = -1,V 2 = -1), and 
Q = 3 (Vi = -1,T>2 = +1). By the symmetry of H AT 
under the interchange a a f-> r a , the sectors Q = 1 and 
Q = 3 are degenerate. Moreover, we observe that the 
ground state belongs to the sector Q — 0. A schematic 
view of the Ashkin- Teller chain is shown in Fig. [3l Note 
that each site contains two spin particles, which means 
that the number N of particles in a chain with M sites 
is N = 2M. 



j-i 



FIG. 3. (Color online) Schematic view of the Ashkin- Teller 
chain. The lattice is composed by two independent spin-1/2 
particles per site j described by Pauli operators {af,Tj"}. 



The model presents an infinite-order quantum critical 
point at j3 = 1 and A = 1. Infinite-order QPTs are 
typically detected by an cxtremum (cither a maximum 
or a minimum) in quantum correlations measures (see, 
e.g., Refs. [1^, [27} for pairwise entanglement and Ref. 0| 
for pairwise quantum discord). However, as shown in 
Ref. (28|, pairwise entanglement is unable to character- 
ize the critical point (/3, A) = (1, 1) in the Ashkin- Teller 
chain. Moreover, it can be shown that pairwise quantum 
discord does not detect such a QPT either. Indeed, tak- 
ing P = 1, the density operator p> ,J '(A) for a pair &j — fj 
is diagonal [28] and exhibits vanishing quantum discord 
for any A. For pairs dj — Sj+i (or fj — Tj+i), the den- 
sity operator ^ ? ' J+1 (A) has off-diagonal terms. Such a 
state displays nonvanishing quantum discord. However, 
no identification (such as an extremum or a cusp) occurs 
at the critical point A = 1 (for any local measurement). 
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FIG. 4. (Color online) GQD associated with the a z eigenbasis 
for a spin quartet in the Ashkin- Teller model for chains up to 
TV = 12 spins. Inset: Derivative of GQD with respect to A. 

On the other hand, if we consider multipartite correla- 
tions, GQD is able to identify the QPT as an cxtrcmum 
at the critical point. However, such identification does 
not occur in the basis that minimizes GQD, which is 
given by the measurement of all spins in the a z eigenba- 
sis. Instead, the infinite-order QPT turns out to be cor- 
rectly characterized if, and only if, local measurements 
are performed in the a x eigenbasis. Remarkably, this is 
exactly the basis of eigenstates of the single spin reduced 
density operators. Therefore, computation of GQD in 
such an eigenbasis can be seen as a generalization of MID 
to the multipartite scenario. 



the <7 2 eigenbasis and in the a x eigenbasis for each par- 
ticle via exact diagonalization of chains up to 16 spin 
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FIG. 6. (Color online) Derivative of GQD associated with 
the <j x eigenbasis for spin quartets, sextets, and octets in the 
Ashkin- Teller model for a chain with TV — 16 spins. 

particles. The results are exhibited in Fig. [4] for mea- 
surements in the a z eigenbasis and in Figs. [5] and [5] for 
measurements in the a x eigenbasis. We can observe that, 
in the a x eigenbasis, the identification of the QPT as an 
extremum (vanishing derivative of the GQD) already oc- 
curs for quartets in lattices with N = 6 spins (see Fig. E|) . 
Moreover, this characterization is kept for sextets and 
octets in larger chains (see Fig. [6J. 

VII. CONCLUSION 
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FIG. 5. (Color online) Derivative of GQD associated with the 
a x eigenbasis for a spin quartet in the Ashkin- Teller model for 
chains up to TV = 12 spins. 

We consider groups of 4 particles (quartets) composed 
by spins &j — <7j+i — Tj — fj+x as well as extensions for 
sextets and octets. For those configurations, we numer- 
ically compute GQD relative to local measurements in 



In summary, we have proposed a measure for multi- 
partite quantum correlations. This measure has been 
obtained by suitably recasting the standard quantum 
discord in terms of relative entropy and local von Neu- 
mann measurements [as given by Eq. (1191) ]. In particu- 
lar, our measure is a systematic extension of the original 
approach for quantum discord as introduced in Ref. Q, 
reducing to it in the particular case of bipartite systems. 
Illustrations of its use have been provided for both the 
Werner-GHZ and the Ashkin- Teller spin chain. Further 
applications of GQD, such as the investigation of mul- 
tipartite correlations in quantum computation and con- 
nections with entanglement (see, e.g., Ref. [HI), are left 
for future research. 
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